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ABSTRACT 


Liquid-filled  holes  used  for  pressure  measurements  of  visco-elastic  liquids 
give  rise  to  systematic  hole  pressure  'errors'.*  Tanner  and  Pipkin  have  presented 
analysis  for  flows  of  a second  order  fluid  in  which  they  derive  a simple  relation 
between  the  first  normal  stress  difference  and  the  hole  pressure  for  flow  situa- 
tions where  Reynolds  numbers  are  very  small.  Implicit  in  the  analysis  is  the 
assumption  that  the  streamlines  are  symmetric  about  the  hole  center  line.  In  this 
paper,  using  a numerical  solution,  we  investigate  the  relationship  between  the 
hole  pressure  and  the  first  normal  stress  difference  for  a range  of  Reynolds 
numbers  where  inertial  effects  are  not  negligible.  The  ratio  of  hole  pressure/first 
normal  stress  difference  is  found  to  vary  from  0.25  to  0.16  as  the  Reynolds  number 
is  varied  from  1 to  10.  Streamline  patterns  are  presented  for  Poiseuille  flow  of 
a second  order  fluid  over  a slot  cut  into  one  wall  of  an  otherwise  straight  channel . 
Various  geometries  are  considered.  The  results  naturally  include  those  for  an  in- 
compressible Newtonian  liquid  at  non-zero  Reynolds  numbers. 

AMS(MOS)  Subject  Classification  70.65  WorkUnitNo.  3-  Applications  of  Mathemat . ic: 
Key  Words:  Visco-elastic  liquids,  Hole  pressure  measurements,  Numerical  Methods 

Permanent  address:  Department  of  Computer  Science,  University  College  of  Swansea,  u.r 

* It  has  been  customary  in  the  past  to  apply  the  term  hole  pressure  'error'  to  the 
quantity  P -P  where  Pij  is  the  pressure  measured  in  the  static  liquid  at  the  end 
of  a hole  in  a channel  wail  and  Pw  is  the  pressure  which  would  be  exerted  on  the  well 
by  liquid  flowing  in  the  channel  if  the  hole  were  not  present.  The  basic  aim  of  thi 
paper  is  to  relate  this  'quantity'  to  material  properties  of  the  liquid,  and  as  such 
the  term  'error' is  misleading  since  t!i«  'quantity'  will  itself  be  subject  to  ei  roi . 

We  choose,  therefore,  to  refer  to  the  'hole  pressure'  defined  as  F -U  . It  will 
be  noted  that  the  sign  of  the  hole  pressure  is  chosen  to  be  posit iveWfor  elastic 
1 i quids . 
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Significance  and  Explanation 

For  same  time  now  the  polymer  processing  industry  has  put  considerable 
effort  into  the  design  of  simple,  reliable  instruments  for  measuring  material 
properties  such  as  viscosity,  elasticity,  etc,  of  polymeric  liquids  or  polymer 
melts. 

This  report  considers  the  theoretical  background  to  one  such  instrument 
which  gives  elasticity  readings  by  measurement  of  pressure  values  in  a liquid 
when  it  flows  through  a straight  channel.  One  wall  of  the  channel  has  a slot 
cut  into  it,  and  pressure  is  measured  at  the  bottom  of  the  slot  and  at  a point 
on  the  channel  wall  immediately  opposite  the  slot.  The  report  looks  in  detail 
at  the  effect  of  inertial  contributions  to  the  flow  and  in  particular  shows 
how  inertial  distortion  modifies  the  relationship  between  the  elasticity  of  the 
liquid  and  the  measured  pressure  values.  Streamline  patterns  are  presented 
showing  how  much  distortion  one  can  expect  for  Reynolds  numbers  up  to  25. 


A COMPUTER  MODEL  OF  HOLE -PRESSURE  MEASUREMENT 
IN  POISEUILLE  FLOW  OF  VISCO-ELASTIC  LIQUIDS 

P.  Townsend 

Introduction 

The  search  for  simple  reliable  instruments  which  measure  material  properties 
of  visco-elastic  liquids  is  a never-ending  one.  In  this  paper  we  use  a numerical 
model  to  investigate  one  flow  situation  in  which  the  use  of  liquid-filled  holes 
or  slots  for  pressure  measurements  of  flowing  visco-elastic  liquids  gives  rise 
to  a direct  measurement  of  the  elasticity  of  the  liquid. 

In  a number  of  very  detailed  and  careful  experiments  carried  out  by  Lodge 
and  his  co-workers  [1,  2 and  3]  it  became  clear  that  errors  are  produced  if 
pressure  measurements  of  visco-elastic  liquids  are  taken  from  transducers  mounted 
at  the  bottom  of  liquid-filled  holes.  Some  work  with  a polvisobutylene  solution 
indicated  an  approximately  linear  relationship  between  the  error  in  the  pressure 
reading  and  the  first  normal  stress  difference.  Lodge  recognized  that  due  to  the 
apparently  systematic  relationship  between  the  'hole  pressure'  and  the  liquid 
elasticity,  a basis  exists  for  an  instrument  which  would  give  direct  readings  of 
elastic  properties  of  a liquid  from  a simple  shear  flow.  These  ideas  are  now 
embodied  in  the  Seiscor-Lodge  Stressmeter  [4] . 

A simple  analysis  due  to  Tanner  and  pipkin  [5]  for  Poiseuille  flow  of  a 
second  order  liquid  over  a slot  showed  that  for  very  small  Reynolds  numbers  where 
inertial  effects  were  negligible,  the  hole  pressure  was  precisely  one  quarter  of 
the  first  normal  stress  difference.  This  result  was  later  confirmed  by  Higashitani 
and  Pritchard  |61  using  a slightly  different  approach.  Kearsley  [7]  derived  a 
similar  result  relating  hole  pressure  and  the  second  normal  stress  difference  for 
slots  placed  alonq  the  main  flow  direction.  The  restrictions  placed  on  these 
theoretical  analyses  were  quite  severe,  however,  and  the  very  simple  results 
obtained  unlikely  to  be  valid  when  fluid  inertia  is  non-negligible.  The  investi- 
gation of  the  effect:;  of  inertia  on  the  relationship  between  the  hole  pressure 
and  the  fir  t normal  stress  difference  in  a visco-elastic  liquid  is  the  aim  of 
this  paper. 
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In  order  to  mathematically  model  a visco~elast ic  liquid,  one  is  faced  with 


1 

the  very  difficult  question  of  which  constitutive  equations  to  vise.  Certainly 
one  would  like  to  be  able  to  model  as  closely  as  possible  as  many  fluid  proper- 
ties as  one  can,  but  often  the  more  sophisticated  fluid  models  qive  rise  to 
partial  differential  equations  which  are  either  too  difficult  or  too  expensive 
to  solve.  In  the  flow  situations  we  shall  consider,  certainly  our  first  concern 
is  for  fluid  elasticity  effects,  and  one  of  the  simplest  models  available  to  us 
is  the  second  order  fluid.  One  has  to  be  careful  in  interpreting  results  obtained 
in  this  case  for,  as  Pipkin  18]  points  out,  the  results  are  valid  only  if  the 
presence  of  second-order  terms  causes  a small  perturbation  of  Newtonian  flow. 

For  flows  which  differ  considerably  from  the  Newtonian  situation  one  may  argue 
that  a second-order  analysis  gives  one  only  preliminary  insight  into  how  Newtonian 
flow  will  be  affected  by  slight  elasticity.  In  the  present  analysis  it  proves  to 
be  the  case  that  the  elastic  solutions  obtained,  except  for  pressure  considera- 
tions, differ  very  little  from  Newtonian  flow. 

The  equations  governing  the  flow  problem,  even  for  the  Newtonian  case,  are 
too  difficult  for  an  analytical  solution,  and  we  must  turn  therefore  to  a numc'rical 
model.  Some  of  the  earliest  work  done  in  this  area  was  that  by  Thom  and  Apelt  |9] 
who  derived  a perturbation  solution  for  small  Reynolds  numbers  for  the  flow  of  a 
Newtonian  liquid.  They  attempted  to  write  down  a simple  relationship  between  the 
hole  pressure  and  the  Reynolds  number.  Some  comparison  of  their  results  is  made 
with  experimental  work  due  to  Ray  l 10] . O'Brien  [11,  l?  and  13]  lias  calculated 
Stokes  flow  of  a Newtonian  liquid  past  slots  of  various  geometries  and  in  situa- 
tions where  the  unperturbed  flow  in  the  main  channel  is  taken  to  be  both  Couette 
flow  and  Poiseuille  flow  and  also  some  combination  of  the  two.  Results  are  pre- 
sented indicating  the  dependence  of  streamline  patterns  on  depth-to-length  slot 
ratios  but  unfortunately  no  reference  is  made  to  hole  pressure  measurement s . 


num  veal  sol"*  ion  of  the  tnl\  Nav  i«»i  -St  okes  oijvi.it  ions  t oi  Newtonian  t'o  i non  i 1 1 o 


t low  past  .»  slot  is  presented  by  Schefenackei  1141,  although  it  is  not  . leal  why 
inertia  t erms  are  ret  aim'll  tot  the  values  of  the  Keynolits  numbet  eons  blot  o.t  ate 
extremely  small.  The  gecmiet  ry  us«>.l  in  this  work  is  derived  tiimi  extension  appli- 
cations amt  i s not  directly  relevant  hole.  Numerical  solutions  foi  high  Keynolits 
number  Newtonian  flow  have  boon  obtained  by  Stevenson  II1'],  who  cons  idols  the  flow 
in  a tube  which  has  a eircumforent ial  wall  cavity  and  calculates  the  stieamlmc 
Again,  howovet  , there  is  no  reference  to  the  hole  prossuio. 

Hole-prossui e calculat  ions  tot  a seoond-oidet  fluid  have  been  made  by  Malkus 
(lt'l,  who  presents  results  f ot  the  rati.'  of  hoi  e-pi  essui  e to  til'd  notmal  sties', 
difference  as  a function  ot  depth-to-widt  h ratio  ot  t he  slot.  both  vViiet  1 1 and 
Poiseuille  flow  are  considered  but  the  analysis  is  limited  to  Stokes  t low. 

In  a recent  paper  by  Crochet  and  bery  (1  *),  con-idei.it  ion  i ••  given  to  the 
flow  .'t  a Maxwell-type  liquid  vn  a u comet  i y similai  to  that  adopted  in  tin-  papei  . 
Only  a limited  numhei  of  tesults  are  pt  osont  rd  since  thv'  problem  is  essential^ 
only  solved  as  one  test  ot  a new  finite  element  technique.  1'hv'  aut  hoi  exj'i  v 
s.-nie  doubt  about  certain  teatuies  ot  theii  solution  due  to  numerical  ditticultie- 
which  they  expel ienoed.  We  shall  make  tuithei  i el  v'i  ence  to  this  woi  k latei. 
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Theoretical  Analysis 


In  this  analysis  we  take  a rectangular  Cartesian  coordinate  system 
(x,y,z)  and  assume  that  we  have  steady  two-dimensional  flow  between  two  infinite 
flat  plates  AH  and  BG  , which  are  positioned  parallel  to  the  0X2  plane  at  a 
distance  h apart  (see  Figure  1) . The  lower  plate  is  assumed  to  have  cut  into 
it  a rectangular  slot  CDEF  parallel  to  the  Z axis,  of  depth  d and  width  b , 
and  of  infinite  length. 

If  we  take  a velocity  vector  v (u,v,0)  then  the  equations  which  govern 


the  notion  of  the  liquid  are 
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where  p is  the  density  of  the  liquid,  P denotes  the  isotropic  pressure  and 

P'  , P'  and  P'  the  relevant  components  of  the  extra  stress  tensor.  The 
xx  xy  yy 

boundary  conditions  which  apply  are  that  the  velocity  components  are  zero  on 
all  plate  surfaces,  and  that  provided  one  is  sufficiently  far  from  the  slot  then 
the  flow  is  undisturbed  Poiseuille  flow.  To  complete  the  specification  of  the 
problem  we  need  to  choose  appropriate  equations  of  state  for  the  liquid.  For 
reasons  given  in  the  introduction  we  confine  attention  to  the  Rivlin-Ericksen 


incompressible  second-order  fluid  given  by 

(1)  . „ _ (2)  ^ . (1)  (l)m 

ik  3 1m  k 


plk  = 2ai  eik'  + 2 ^ e;:#  + 4 a’  e;:  e* 
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where  a,,  a_  and  a,  are  material  constants  and  e.:"  is  the  jth  rate  of 
12  3 ik 

strain  tensor  introducted  by  Oldroyd  [18]^.  This  fluid  exhibits  both  first  and 
second  normal  stress  differences,  governed  by  non-zero  values  of  the  parameters 
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where  a!,^  is  the  jth  Rivlir.  - Tricksen  tensor, 
ik 
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but  is  restricted  to  a constant  viscosity 


a,  and  , but  is  restricted  to  a constant  viscosity  . 

The  throe  relevant  extra  stress  components  are  then  given  by 
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We  now  take  a number  of  steps  to  simplify  the  governing  equations.  Kitst 
we  introduce  the  following  set  of  non-dimensional  variables 


v*  - I 


+ - V 

. u*  - - , V*  - - 


t.  ' 1.  V ' V 

\ (ttl 

r 2 t ( 

b*  - — , a*  - — , a*  - — 
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where  l.  is  some  charact  er  ist  ic  length  and  l'  is  a chat  act  er  i st  io  velocity.  Koi 
simplicity  we  choose  immediately  to  drop  the  * notation  although  it  is  still 
imp l ied . 

Secondly  we  introduce  a stream  function  v(x,yl  such  that 

u - -p*-  and  v - - v*  • (v>3 

,ly  ,'x 

With  this  definition  the  equation  of  continuity  ('1  is  satisfied  identically. 

finally  we  substitute  the  expressions  (M-(71  fo»  the  stress  component  s 
into  egu.it  ions  (11  and  (.'1  and  eliminate  the  pressure.  The  final  equations  which 
we  then  obtain  are 


_1_ 

Re 


(10) 


Simplified  Analysis  of  Tanner  and  Pipkin  [5] 

If  one  assumes  that  the  flow  is  so  slow  that  one  may  neglect  inertia,  then 
the  disturbance  due  to  the  slot  is  symmetrical  about  the  slot  centerline,  and 
therefore  the  pressure  p in  a Newtonian  liquid  is  constant  along  this  line. 

If  one  further  assumes  that  the  slot  is  sufficiently  deep  that  there  is  negli- 
gible motion  at  the  bottom  of  it  then  this  constant  value  is  the  pressure  P() 
at  the  bottom  of  the  slot.  Under  these  restrictions  Tanner  and  Pipkin  show  that 

for  a second  order  fluid  the  stress  along  the  centerline  is  given  by 

„ „ x n (1)  . , (1)  (l)m  1,.  , (1)  (l)m..  ,, 
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+ a„[2e.,  - ( — — - 3 trace  (e.  e.  ))<5..J 

3 ik  DT  im  k lk 

If  one  now  further  assumes  that  the  slot  in  the  bottom  plate  is  so  narrow 

(h/b  large)  that  the  flow  near  the  top  plate  is  negligibly  disturbed  then  it  is 

fairly  straightforward  to  show  that  the  centerline  thrust  P on  the  top  plate 

w 


is  given  by 


p - pH  - i «.(  f >2 

w H 2 2 3y 


Since  the  first  normal  stress  difference  N^  is  given  by 

N ,=  P - P = -2a_ ( )2  (17) 

l xx  yy  2 3y 

then  from  (16)  and  (17)  we  have  the  simple  relationship  between  the  first  normal 

stress  difference  and  the  hole  pressure  P - P„  namely 

w H 

P - P = 7 N . (18) 

w H 4 1 

Ttiis  result  suggests  that  measurement  of  the  hole  pressure  is  likely  to 


give  a direct  measure  of  the  elasticity  of  a liquid.  Ttiis  analysis  tias  been 
extended  by  Kearsley  (7)  who  was  .able  to  show  that  for  rectilinear  flows  along 
a slot  the  hole  pressure  is  one  half  the  second  normal  stress  difference. 
Higashitani  and  Pritchard  [6)  using  a somewhat  different  approach  have  confirmed 
the  Tanner-Pipkin  result  although  here  again  the  methods  used  depend  crucially  on 


the  flow  patterns  being  symmetric  about  the  centerline  of  the  slot.  This 


is  a severe  restriction.  In  practical  instruments  one  has  a finite  Reynolds 
number  and  one  must  expect  inertial  effects  to  produce  asymmetrical  stream- 
lines within  the  slot.  To  investigate  fully  the  relationship  between  hole- 
pressure  and  the  first  normal  stress  difference  one  must  retain  the  inertia 
terms  in  Equation  (11)  and  derive  solutions  for  some  finite  positive  Reynolds 
number.  In  this  case,  however,  analytical  solutions  are  out  of  the  question 
and  one  must  turn  to  numerical  techniques. 


Numerical  Solution 


In  order  to  solve  the  equations  numerically  it  is  convenient  to  introduce 
an  intermediate  variable  f,  as  follows.  We  write  Equation  (10)  in  the  form 


_L_  v2r_iiii_Mi£=r 
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and 
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£(x,y)  represents  a deviation  from  Newtonian  behaviour.  If  is  zero  then, 

from  (20)  , E,  is  identically  zero  and  (11)  and  (19)  reduce  to  the  Navier-stokes 
equations . 

To  solve  these  equations  together  with  (11)  , we  use  finite  difference  tech- 


(21) 


To  discretize  the  Laplacian  operators  in  (11)  and  (19)  we  use  standard  five  point 
difference  formulae 


; mesh  (x^,y^) 

(see  Figure  2)  defined  by 

x.  = i Ax  , 
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For  the  non-linear  terms  in  (19)  it  is  necessary  to  *-ake  special  measures  to  insure 
diagonal  dominance  of  the  difference  equations.  Otherwise  it  has  been  found  that 
the  iterative  scheme  to  be  applied  to  solve  these  equations  will  not  converge  for 
Reynolds  numbers  in  excess  of  about  unity.  One  technique  used  by  Schafenacher  (14) 
is  to  introduce  additional  terms  into  the  difference  approximation  for  the  x 
derivative  of  f,  , based  on  the  previous  iterative  step.  Another  more  popular 
technique  is  to  make  use  of  so-called  'upwind'  differencing.  In  this  case,  one- 
sided difference  approximations  to  the  first  derivatives  of  l are  used,  the  parti- 
cular approximation  i.e.  forward  or  backward,  being  chosen  according  to  the  signs 


-10- 


of  the  corresponding  derivatives  of  1J1  , to  ensure  diagonal  dominance.  In  the 


work  considered  here',  the  latter  technique  has  proved  to  be  a somewhat  more 
satisfactory  method  although  one  has  to  accept  a loss  of  accuracy  due  to  the  first 
order  nature  of  the  one-sided  difference  approximations.  Greater  accuracy  can  be 
recovered  however,  if  required,  by  applying  a difference  correction  procedure 
(See,  for  example,  Dennis  and  Cheng  (191). 

Discretization  of  Equation  (20)  presents  similar  diagonal  dominance  problems 
and  we  must  again  use  one-sided  differences  for  the  first  derivatives  of  l,  . 

For  boundary  conditions  we  have  that  i)'  and  its  first  derivatives  are  zero 
on  solid  boundaries,  and  for  the  inlet  and  outlet  conditions  we  impose  Poiseuille 
flow  profiles.  As  an  alternative  condition  at  the  outlet,  <Ji  may  be  assumed  to 
be  independent  of  x , althouqh  this  slows  down  the  convergence  of  the  iterative 
procedures.  For  the  range  of  Reynolds  numbers  considered  here,  it  proves  adequate 
to  impose  both  inlet  and  outlet  profiles. 

The  three  difference  equations  may  now  be  solved  using  successive  over- 
relaxation, although  because  of  the  coupled  nature  of  the  system  it  is  necessary 
to  set  up  an  inner/outer  iteration  procedure.  Normally  to  start  such  a procedure 
one  makes  initial  guesses  ij>  , t,  and  f,  for  the  interior  of  the  finite 

difference  mesh,  and  then  using  an  appropriate  Taylor  Series  formula  one  constructs 
from  boundary  values  for  f.  ^ . (See,  for  example,  Greenspan  and  Schultz 


120]).  From  t, 

, „ • , (1) 
i/'  . Using  i/i 

The  sequence 


^ one  may  then  solve  (11)  to  give  a new  approximation  for 


and  f.  ^ one  then  proceeds  to  solve  (19)  for  t ^ , and  so  on. 


(23) 
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of  approximations  (outer  iterates)  to  each  of  the  three  variable.',  is  calculated 
until  the  nth  and  (n+l)st  outer  iterate  agree  to  some  predetermined  tolerance. 


In  some  cases  convergence  problems  are  experienced  and  it  is  found  necessary  to 
apply  smoothing  between  each  outer  iterate.  When,  for  example,  .1  new  approxima- 
tion is  computed,  then  a weighted  mean  ij/  of  this  value  and  the 


previous  value  given  by 

7 (n+1)  (n)  . , (ntl)  , . .. 

4'  = p„  4'  + (1  - p ) 4’  , 0 < P v 1 (^.4) 

2 s 

is  actually  used  for  4'  for  the  next  stage  of  the  outer  iteration  procedure.  A 
similar  smoothing  is  applied  to  r.  and  to  f.  . 

Greenspan  and  Shultz  (20)  have  suggested  that  it  may  be  possible  to  avoid 
this  inner/outer  iteration  procedure.  They  solve  a number  of  problems  for  the 
Navier  Stokes  equations  in  which  successive  over-relaxation  is  applied  simultane- 
ously to  the  vorticity  and  stream  function  equations.  They  found  that  this  was  a 
much  faster  procedure  which  lias  the  added  advantage  of  not  requiring  smoothing 
parameters . However,  attempts  to  implement  these  ideas  lei  the  problems  con- 
sidered here  proved  only  to  be  partly  successful.  Simultaneous  iteration  ot  all 
three  equations  (11),  (19)  and  (20)  did  not  produce  a convergent  solut  ion,  and 
even  if  only  the  first  two  equations  were  iterated  tixjethei  t hen  a convergent 
solution  is  only  obtained  if  a certain  degree  of  smoothing  is  applied.  in  this 
case,  a smoothing  correction  is  applied  at  the  end  of  each  complete  sweep  through 
the  mesh  of  the  41  and  iterations.  This  method  did,  however,  prove  to  be 

much  faster  than  an  inner/outer  iteration  scheme.  4 and  ate  now  treated  as 

a pair  and  combined  with  f in  a normal  inner/out ci  it erat  ion  as  follows 


(0)  (0) 


.(0) 


(1)  (11 


' L 


(.'1  (.') 
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In  spite  of  the  efforts  set  out  above,  a convergent  solution  proved  difficult 
to  obtain  for  Reynolds  numbers  in  excess  of  about  SO.  With  even  greater 
smoothing  and  the  use  of  a finer  finite  difference  mesh  it  is  likely  that  the 
range  of  Reynolds  numbers  could  be  extended.  Kor  our  purposes,  we  are  most 
interested  in  somewhat  slower  flows,  and  therefore  no  attempt  was  made  to 
extend  the  analysis  to  very  high  Reynolds  numbers. 


r 


■ ■ 


The  Hole  Pressure 

It  is  difficult  to  define  precisely  what  we  mean  by  the  hole  pressure. 

Essentially  if  a pressure  transducer  is  flush  mounted  in  the  bottom  wall  of  the 

slot  and  gives  a reading  P , and  a second  transducer  mounted  directly  opposite 

H 

in  the  top  plate  gives  a readinq  P then  the  difference  P -P„  is  the  hole 

w w H 

pressure . 

However  difficulty  arises  as  to  how  to  simulate  the  action  of  the  transducer 
diaphragm.  It  was  decided  to  compute  an  average  force  exerted  by  the  fluid  on 
seme  portion  of  the  bottom  wall  of  the  slot  together  with  the  equivalent  average 
force  exerted  on  a similar  portion  of  the  top  plate.  The  range  of  integration 
chosen  in  the  averaging  was  the  complete  bottom  wall  of  the  slot  except  for  a 
portion  .\x  at  each  end.  This  avoids  diff iculties  at  the  corners  C and  F . 
The  hole  pressure  A P is  then  given  by 


A P » /(P  - P )dx/(b-2Ax) 

^top  plate  ^slot 

whetc  p i'.  t lie*  normal  component  oif  the  stress  tensor  given  by 

yy 

p - -p  t p' 

yy  yy 

The  difference  in  the  values  of  the  fluid  pressure  P is  given  by 


(26) 


(27) 


P - p 

top  plate  slot 


where  from  (2) , (b)  and  (7) 
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Both  inteqrat ions  indicated  above  are  carried  out  numerically  using  a simple 
Simpson's  rule  algorithm. 

In  order  to  compare  our  result  with  the  Tanner-Pipkin  result,  a final  ratio 

K , corrected  for  inertial  effects,  is  computed,  where 

K = (AP  - AP  )/N  (291 

second  order  fluid  Newtonian  ' 1 K ’ 

and  where  is  the  first  normal  stress  difference  defined  as 

N = P - P . (30) 

1 xx  yy 

To  compute  we  take  average  values  of  the  stress  components  over  the  region 

of  the  top  plate  opposite  the  slot. 

For  the  different  slot  geometries  considered,  different  finite  difference 

meshes  were  constructed.  In  all  cases,  however,  approximately  20  qrid  cells  wort 

chosen  to  represent  the  y variation  of  solution  in  the  main  portions  ot  the 

channel  away  from  the  slot,  and  likewise  approximately  20  grid  cells  wet e chosen 

to  represent  the  x variation  within  the  slot.  The  over-all  total  of  grid  point 

was  dependent  on  the  particular  geometry  under  considcrat  ion . .Ml  eomput at  ion  ... 

carried  out  on  a ItNIVAC  1110  computer. 
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Result  s 


Streamlines  for  various  flow  conditions  are  shown  in  Figures  3-8.  It  will 
be  noted  that  in  ail  cases  a secondary  flow,  very  much  weaker  than  the  flow  in 
the  mainstream,  is  set  up  in  the  slot. 

For  creeping  flow  it  is  well  known  that  the  velocity  field  of  a Newtonian 
liquid  is  identical  with  that  of  a second  order  fluid.  See  for  example  diesekus 
[211,  Tanner  [22].  It  is  of  interest  to  investigate  just  how  much  the  twit  flow 
fields  differ  when  inertial  effects  are  included.  Since  the  first  normal  stress 
difference  acts  like  a tension  along  the  streamlines  that  tends  to  pull  the  liquid 
out  of  the  slot  one  would  expect  that  part  of  the  Newtonian  flow  field  which  dips 
into  the  slot  to  straighten  for  a second  order  fluid.  Certainly  for  the  range 
of  Reynolds  numbers  considered  here  any  changes  due  to  elasticity  proved  to  be 
in  the  expected  direction  but  these  changes  are  extremely  small  in  magnitude. 

This  observation  was  also  made  by  Crochet  and  Bezy  [17] . In  Figure  3 we  have 
plotted  streamlines  for  a square  slot  of  depth  equal  to  the  main  channel  height. 
The  full  lines  are  for  a Newtonian  liquid  at  a Reynolds  number  of  25  and  the 
dashed  lines  indicate  changes  due  to  elasticity  when  a*  = -.1.  It  can  be  seen 

4- 

lust  how  little  the  flow  field  is  modified  by  elasticity  even  thouqh  inertial 
effects  are  suf f iciently  large  to  cause  the  overall  flow  pattern  to  deviate  quite 
considerably  from  the  pattern,  symmetric  about  the  centerline,  which  one  obtains 
for  creeping  flow. 

In  Figures  4 and  5 we  can  see  how  the  gradual  influence  of  inertia  is  felt 
as  the  Reynolds  number  increases.  At  a very  small  Reynolds  numbei , the  stream- 
lines across  the  whole  channel  are  modified  by  the  presence  of  t lie  slot,  those 
close  to  the  slot  being  drawn  quite  deeply  into  the  slot.  As  the  Reynolds  numbet 
increases  the  streamlines  straighten  until  only  those  quite  close  to  the  slot  are 
affected.  Another  noticable  effect  of  increasing  the  speed  of  the  mainstieam  is 
I o strengthen  t he  secondary  flow  in  the  slot  although  it  always  remains  several 
orders  of  magnitude  weaker  than  the  main  flow. 


t 
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Figure  4 


Streamline  projections  for  a square  slot  with  h/b  = 1 


and  for  a second  order  fluid  (a*  = -.1) 


Full  lines  - Re  = 1.  Broken  lines  - Re  * 10. 
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In  Figures  6-8  we  have  considered  how  the  flow  field  is  affected  by  changes 
in  the  channel/slot  geometry.  Two  factors  are  relevant  here,  namely  the  slot 
depth  and  the  slot  width.  Figures  6 and  7 show  streamlines  plotted  for  slots 
of  half  and  double  depth  respectively.  A comparison  of  Figures  4,  6 and  7 shows 
that  the  main  flow  is  relatively  unaffected  by  slot  depth  and  that  in  all  cases 
the  fluid  is  virtually  stationary  at  the  bottom  of  the  slot.  One  might  expect 
therefore  that  hole  pressures  would  be  relatively  insensitive  to  slot  depth  and 
we  shall  investigate  this  point  later.  One  consequence  of  changing  the  slot  depth, 
however,  is  a considerable  change  in  the  secondary  flow,  and,  in  particular,  for 
slots  of  double  depth,  one  sees  the  onset  of  a second  separate  circulatory  motion 
in  the  slot. 

Figure  8 shows  streamlines  for  a slot  of  half  width.  Intuitively  one  might 
expect  inertial  effects  to  be  reduced  for  a narrower  slot  and  the  results  seem  to 
confirm  this,  for  the  streamlines  in  the  main  channel  are  straighter  compared 
with  a square  slot,  the  overall  pattern  is  more  symmetric  about  the  slot  center- 
line,  and  very  little  change  is  seen  when  the  Reynolds  number  is  increased  from 
10  to  25. 

If  we  turn  our  attention  now  to  the  hole  pressure.  Figures  9 and  10  show 
plots  of  the  hole  pressure  as  a function  of  distance  across  the  slot.  Further 
tabulated  values  for  representative  points  across  the  slot  are  calculated  in  Table 
1.  It  can  be  seen  from  the  figures  that  for  the  Newtonian  case  (dashed  lines) , 
a positive  value  in  the  first  part  of  the  slot  is  to  some  extent  cancelled  by  a 
negative  value  in  the  remainder  of  the  slot  so  that  the  net  integrated  hole  pres- 
sure recorded  by  two  transducers  will  be  very  small.  This  is  not  the  case  for 
the  elastic  liquid  (full  lines)  where  the  curves  remain  positive  throughout  the 
slot.  One  would  expect,  therefore,  the  hole  pressure  contribution  due  to  elas- 
ticity to  be  at  least  an  order  of  magnitude  greater  than  that  due  to  inertial 
effects  at  the  Reynolds  numbers  considered  here. 
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Hole  pressure  plotted  as  a function  of  dist; 
across  the  slot  for  a narrow  slot  with  h/b 


Although  the  parameter  does  not  appear  in  the  equations  for  the  stream- 

lines, it  does  affect  pressure  calculations.  Throughout  this  work  a variation  ot 
n.j  within  an  acceptable  range  of  values  proved  to  have  a virtually  negligible 
effect  on  the  results. 

In  Figure  11  we  have  plotted  the  ratio  R as  a function  of  Reynolds  number. 

It  is  clear  that  the  assumptions  made  in  the  Tanner-Pipkin  theory  breakdown 
fairly  rapidly  as  the  Reynolds  number  increases  above  unity  although  there  is  still 
a very  close  relationship  between  the  hole  pressure  and  the  first  normal  stress 
difference.  Also  plotted  in  Figure  11  are  some  of  the  results  of  Crochet  and 
Brezy.  It  should  be  said  that  it.  is  not  clear  from  their  work  whether  they  have 
applied  an  inertial  correction  as  we  have,  by  subtracting  off  the  Newtonian  contri- 
bution to  the  hole-pressure.  (See  Equation  (29)).  This  correction  would  brinq 
their  results  a little  closer  to  those  computed  here  although  some  difference  is 
still  apparent. 

Finally  we  return  to  the  question  of  geometrical  effects  in  hole  pressure 
measurement.  In  Table  2 the  ratio  R is  tabulated  for  four  different  slot 
geometries.  As  suggested  earlier,  it  would  seem  that  the  relationship'  between 
hole  pressure  error  and  the  first  normal  stress  difference  is  relatively  insensi- 
tive to  changes  in  the  depth  of  the  slot.  The  width  of  the  slot  is  a more  im- 
portant factor  however.  A narrower  slot  leads  to  less  inertial  distortion  of  the 
flow  field  and,  as  a result,  a value  of  K closer  to  that  predicted  by  the 


Tanner-Pipkin  theory. 


Half  depth 


Double  depth 


Half  width 


Table  2 


The  ratio  R tablulated  as  a function  of 


Cone  \ us 1 •. wu; 

The  present  invest  iqat  ion  into  the  determination  of  the  first  normal  stress 
difference  of  a visco-elastic  liquid  by  the  measurement  of  hole-pressures  indi- 
cates that  the  probltm  is  somewhat  more  complex  than  that  modelled  by  the 
Tanner /Pipkin  theory.  Although  elasticity  is  found  to  have  very  little  influ- 
ence on  streamline  patterns,  even  at  quite  appreciable  Reynolds  numbers,  inertial 
distortion  of  the  flow  field  results  in  considerable  change  in  the  relationship 


between  the  hole-pressure  and  the  first  normal  stress  difference  when  Re  N 1 . 
However,  our  results  confirm  the  Tanner/Pipkin  theory  for  Reynolds  numbers 
below  unity. 

Slot  geometry  is  not  found  to  be  an  important  factor  as  far  as  slot  depth 
is  concerned,  but  narrower  slots  are  found  to  reduce  inertial  distortion  of  the 


The  results  confirm  that  elastic  and  inertial  contributions  to  the  hole 

pressure  are  additive  in  the  sense  that  AT,  , . - AP  . is  propor- 

elastic  Newtonian 

tional  to  N ^ for  given  Re  and  h/b  . 
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Abstract  (continued) 


s 

Reynolds  numbers  where  inertial  effects  are  not  negligible.  The  ratio  of  hole 
pressure/first  normal  stress  difference  is  found  to  vary  from  0.25  to  0.16  as  the 
Reynolds  number  is  varied  from  1 to  10.  Streamline  patterns  are  presented  for 
Poiseuille  flow  of  a second  order  fluid  over  a slot  cut  into  one  wall  of  an  other- 
wise straight  channel.  Various  geometries  are  considered.  The  results  naturally 
include  those  for  an  incompressible  Newtonian  liquid  at  non-zero  Reynolds  numbers. 


It  has  been  customary  in  the  past  to  apply  the  term  hole  pressure  'error'  to  the 
quantity  p^-p  where  is  the  pressure  measured  in  the  static  liquid  at  the  end 

of  a hole  in  a channel  wall  and  P^  is  the  pressure  which  would  be  exerted  on  the 

wall  by  liquid  flowing  in  the  channel  if  the  hole  were  not  present.  The  basic  aim 
of  this  paper  is  to  relate  this  'quantity'  to  material  properties  of  the  liquid, 
and  as  such  the  term  'error'  is  misleading  since  the  'quantity'  will  itself  be  subject 
to  error.  We  choose,  therefore,  to  refer  to  the  'hole  pressure'  defined  as  P -P  . 

It  will  be  noted  that  the  sign  of  the  hole  pressure  is  chosen  to  be  positive  for 
elastic  liquids. 


